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Abstract. We present a closest separable state to cluster states. We start by 
considering linear cluster chains and extend our method to cluster states that can be 
used as a universal resource in quantum computation. We reproduce known results for 
pure cluster states and show how our method can be used in quantifying entanglement 
in noisy cluster states. Operational meaning is given to our method that clearly 
demonstrates how these closest separable states can be constructed from two-qubit 
clusters in the case of pure states. We also discuss the issue of finding the critical 
temperature at which the cluster state becomes only classically correlated and the 
importance of this temperature to our method. 
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1. Introduction 

Entanglement plays an important role in modern physics and is the subject of intense 
research both for its implications in our fundamental understanding of nature pQ as well 
as its practical usefulness in quantum information processing [2] and [3]. Different 
models of quantum computing have been devised to harness this usefulness. The 
most prominent ones include the circuit model [I], topological computing [5] and 
measurement-based computing [6]. 

Measurement-based quantum computation (or one-way quantum computation) 
differs from the circuit model by using an initially prepared highly entangled resource 
state. The algorithm proceeds as a set of local measurements on this state. The basis 
and order of these measurements characterize the algorithm and outcomes of previous 
measurements are fed forward to make the computation deterministic. 

The general resource for this computation model are graph states [TJ. In this paper 
we will concentrate on particular graph states, namely cluster states [H]. Unlike graph 
states, cluster states have regular square structure. One of the advantages of one-way 
quantum computation is that all the entanglement used in the algorithm is present 
right at the beginning. The subsequent measurements only consume it. This makes 
it easier to identify the role of entanglement in the information processing compared 
to the traditional circuit model. Recently it has also been shown both theoretically 
and experimentally in [9] that the quantum computational power of cluster states 
can be extended by replacing some projective measurements with generalized quantum 
measurements. This makes cluster state quantum computation a very attractive and 
promising prospect when it comes to constructing a quantum computer. 

Due to its clear role as a resource in one-way quantum computation it is highly 
desirable to know the entanglement properties of cluster states. Perhaps one of the 
most important problems in studying the properties of quantum states, and resource 
states for quantum computation in particular, is the quantification of entanglement 
itself. The entanglement scaling of pure d- dimensional graph states has been shown by 
Markham et al to be N/2 [10]. However any realistic quantum computer will operate at 
finite temperatures which makes it necessary to be able to characterize entanglement of 
mixed states. Quantifying entanglement is a notoriously difficult task [UJ and the exact 
scaling of entanglement in thermal cluster states is still unknown though some advances 
have been made by investigating the localization of entanglement in noisy cluster states 

H2|. 

There are many entanglement measures that quantify the scaling of entanglement 
in quantum states. In this work we use relative entropy of entanglement (see [TTJ and 
[18J) to quantify the scaling. In general there are not many multipartite quantum 
states for which the relative entropy of entnaglement can be calculated analytically. 
D-dimensional symmetric states are one of the few examples [13] . |14j . It is unknown 
how to find the closest separable state to a general entangled state and the problem of 
closed form of relative entropy still remains an open question. Though some progress 
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has been made in this direction ( [15] and [TH]) it is impossible to do this even for a 
general two-qubit state. 

We present a method of constructing the closest separable state to a thermal cluster 
and show how it can reproduce the results of [10J for pure cluster states. We also show 
how bound entanglement complicates the procedure of finding the closest separable state 
for noisy cluster states. 



2. Closest separable state for pure cluster states 

The entanglement of cluster states is known to scale linearly with the size of the system 
as y. This result was first obtained by Markham et al in [TU| . The authors used a 
technique where they derived upper and lower bounds on the entanglement and then 
they showed that these two bounds are equal. 

In this section we verify this result using geometrical ideas. To do this we employ 
the relative entropy of entanglement ([IT]. [IS]) E(a) defined as 

E(a) = mm S(a\\p) (1) 

where a is the entangled state. The minimization is taken over the set D of all separable 
states p and S(a\\p) = Tr[oTog cr — crlogp] is the relative entropy. For pure states this 
simplifies to S(a\\p) = — Tr[alogp]. Since we are considering only qubits we take the 
logarithm to base 2. State that minimizes the relative entropy in (JTJ) will be designated 
by p* and will be called the closest separable state. In order to prove that this state p* 
achieves the minimum of relative entropy we employ the same approach as [19] . Consider 
a small region around the closest separable state (1 — x)p* + xr for a small parameter x. 
All we have to show now is that for a general separable state r the gradient of relative 
entropy is non-negative in this region. Due to convexity of the set of separable states 
this will also mean that the minimum of relative entropy is global. 

A cluster state can be thought of as a state whose qubits are first prepared 
in the +1 eigenstate of a x and then control-phase gates CZ are applied between 
nearest neighbours. For example the two-qubit cluster state is \ip 2 ) = CZ\+)\+) = 
|[|00) + 1 1 } + 1 10} — |11)]. It can be easily seen that the closest separable state p\ is 

/ 1 1 \ 

P*- 1 - 1 1 ° ° (2) 
P2 " 4 1 -1 {2) 

K -1 1 i 

by calculating that E(a 2 \\p2) = 1 where a 2 = \ip2) (^2 1- 

So far all the states have been expressed in the computational basis {|0), |1)}. 
It is however more useful to write the states in what we call the "mixed" basis 
{|0), |1), |+), |— )}. So the state vector and the closest separable state for the case of 2 
qubits now become 

hk> = ^[|0+) + |l->] (3) 
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Pi =^[|0+)(0+| + |l-)(l-|] (4) 
Similarly for 4 qubits the expressions take the following form 

|V> 4 ) = i[|0 + 0+) + |0 - 1-) + |1 - 0+) + |1 + 1-)] (5) 

Pi = -[|0 + 0+)(0 + 0+| + |0-l-)(0-l-| (6) 

+ |1 - 0+)(l - + | + |1 + + 1 - |] 

From the form of p* 2 and p\ we can see that the closest separable state is obtainable from 
the state vector by writing the density matrix in the "mixed" basis and then keeping 
only the diagonal terms. This is why it is more instructive to work in the "mixed" basis, 
one can immediately see the form of the closest separable state from the state vector. 
Of course in the computational basis not all off-diagonal terms are zero. 

Important thing to notice is the number of terms in the state vector of a cluster state 
is different in the computational basis and in the "mixed" basis. In the computational 
basis all the coefficients are non-zero, they are either 1 or —1. So for iV-qubit cluster state 
there are 2^ basis coefficients. On the other hand in the "mixed" basis the number of 
coefficients drops to 2 iV//2 . So when it comes to normalization coefficients for the density 
matrices we have 





Vector \iI>n) 


Cluster a N 


Separable p* N 


Computational basis 


2-JV/2 




2-n 


"Mixed" basis 


2~N/4 


2-N/2 


2-N/2 



From the above it is clear to see that the closest separable matrix p* N is of rank 2 N I 2 . 
Another useful fact that will be used later is that when the matrices are expressed in the 
"mixed" basis it is easy to see that they commute with the cluster state [p* N , a N ] = 0. 

Before we start the proof for the closest separable state a small note about our 
notation is in place. The proof contains places where in stead of just writing or 
p* N it is more instructive to write out the state explicitly, for instance for 2 qubits 
a = |(|0+) + |1— ))*((0 + | + (l — |). This is of course not possible in the case of general 
N. Therefore we use the following notation. 

Pure vector \ip N ) = 2- N/ \\...) + ...) 

Pure matrix a N = 2 _Af/2 (|...) + ...) * ((...| + ...) 

Closest sep. state p* N = 2- N/2 {\...)(...\ + ...) 

From the normalization factors it is clear that the states are in the "mixed" basis. This 
also implies that each round bracket contains 2 N I 2 terms. 

As mentioned above, to prove that states of the form p* 2 and p\ are really the closest 
separable states we need to show that the gradient of relative entropy in a small region 
around p* N is non-negative 

\im^-S(a\\(l-x)p* N + xr)>0 (7) 
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Substituting f(x,r) = S(a\\(l — x)p* N + xr) and realizing that for a positive operator 
A we have the expression log A = J °° ^j^-ppp condition ([7]) becomes 

|£(0, r) = Tr[a jf" \p* N + t)-\p* N - r)(p* N + t^dt] > (8) 
Starting with the first integral in Eq.® 

POO 

= Tr[a / (p^ + ty'pM + tr'dt] 
Jo 

POO 

= Ti[a (p* N + t)- 2 dtp* N } (9) 
Jo 

Evaluating the integral in Eq.(j9]): 

) POO 

( P * N + t)- 2 dt= / (2-*/ 2 + ty 2 dt(\. ..)(...| + ...) 



(2- N ' 2 + t)- 2 dt*2 N ' 2 p* N 
= 2 N ' 2 * 2 N / 2 p* N 

= 2 N P*n (10) 
Substituting this back into Eq.Q: 

Tr[a * 2*^] = 2 N [a * 2^(|...)(...| + ...) * (|...)(...| + ...)] (11) 

= Tr[a(|...)(...| + ...)] 

= 2- N ' 2 * 2 N ' 2 = 1 

Now turning our attention to the second integral in Eq.fjHJ). Here we will use the fact 
that a and p* N commute with each other. 

POO 

= Tr[/ (p^ + t^aip^ + ty'dtr) 
Jo 

/oo 
(p* N + t)- 2 dtar] 

= 2 N Tr\p* N or] 

= 2 Af Tr[2- Ar / 2 2" iV / 2 (|...)(...| + ...) * (|...) + ...) * ((...| + ...)r] 

= '&[(|...><...| + ...)*(|...) + .. .)*((...|+.. .)r] 
= Tr[(|...) + ...)*((...| + ...)r] 

= 2 N / 2 [ar] (12) 
So finally the gradient of the relative entropy in Eq.fjHJ) can be written as 

£(0,r) = l-2 Af / 2 Tr[ar] (13) 

In order for inequality ({TBI to be satisfied the trace needs to scale at most as 
2~ N ' 2 . This is in fact the case as can be proved by induction. First we will consider 
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Figure 1. Plot of 2Tr[<72T2] showing that the maximum of the trace is 1 and therefore 
the gradient of relative entropy is non-negative. Here a and b are parameters from r 
representing the different amplitudes for the two subsystems. So r = \a(3){a(3\ where 
| a) = a|0) + y/l - |a| 2 |l) and similarly \0) = 6|0> + y/1 - \b\ 2 \l). 



t = \a/3)(a/3\. The inequality has been verified numerically for 2, 4 and 6 qubits. In 
the case of 2 qubits graphical solution has been obtained as well in Fig.([T]). 

N = 2 Tr[(7 2 r 2 ] <\ = ^ (14) 

N = 4 Tr[a 4 r 4 ] < ^ = T 2 

N = 6 Tr[a 6 r 6 ] < i = 2~ 3 

o 

Assuming that the above holds true for the case of N = k, namely Tr[<7fcTfc] < 2~ fc / 2 . 
Now we need to prove using the above assumption that the case N = k + 2 holds true 

Tr[a k+2 r k+2 ] < 2- k ' 2 - 1 (15) 

Starting with the trace 

Tr[cr fc+2 r fe+2 ] = (^ k+2 \r k+2 \ij; k+2 ) (16) 
= (<iP k \(ip 2 \C&T k ®T 2 CZ\iP k )\< l p 2 ) 

where we have used the fact that l^fc+2) = CZ\ip k )\ip 2 ) and r k+2 = r k ®t 2 . The control- 
phase gate acts on qubits k and k + 1. Also the separable state r does not commute 
with the controlled phase gate so their commutator is non-zero \r k ®T 2 ,CZ\ = A so we 
can substitute r k <g> t 2 CZ = A + CZr k ® r 2 into Eq. flTB]) . 

Tr[a fe+2 r fc+2 ] = (^| (^|CZU|^) |^ 2 ) (17) 

+ MTkM * (^2|T 2 |^ 2 ) 
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Taking a closer look at the product CZ^A and using the above commutation relation 
we get CZ^A = CZ^t^ <g> t 2 CZ — <g> r 2 . Taking the trace of both sides 

Tr[CZU] = Tr[CZV fc ® r 2 CZ] - Tr[r fc ® r 2 ] = (18) 

Using the assumption for N = k this implies 

Tr[<7 fc+2 T fe+2 ] = (^|r fe |^) * (H^M < 2~ k/2 - x (19) 

Therefore (df/dx)(0, \a(3 . . .){a(3 . . . |) > 0. Since any separable state can be written in 
the form p = '^2 li p l \ct L (5 l . . .){a l /3 l ... | we have 

§£(0, P) = E^ ' 1^ • • • • • I) ^ ( 2 °) 

i 

Therefore p* N is the closest separable state to a^. Now we can compute the relative 
entropy between the cluster state and p* N . 

E(a N ) = S(a N \\p* N ) (21) 

= -Ti[a N \ogp* N ] 

= - (V^llogp^Tv) 

= -\og(ip N \p* N \ip N ) 

= - \o g 2- N ^2- N / 2 ((...\ + ...) * (|...)(...| + ... + |...) * (|...> + ... + |...) 

= - log 2- N 2 N ' 2 

= -log2-^ 2 = ^ (22) 

This concludes the proof that our closest separable state gives linear scaling of 
entanglement and reproduces the results of [TO] . 

We can look at the form of the closest separable state from a more operational 
point of view. Imagine we are asked to prepare the closest separable state to an iV- 
qubit one-dimensional cluster state. The only ingredients needed to construct this state 
are the 2-qubit closest separable state p* 2 and control-phase gates CZ. The state p* N can 
be written 

p* N = U * {pi ® p* 2 ® . . . <g> p* 2 ) *U ] (23) 

V v ' 

N/2 

where the unitary applied has the form 

U = I®CZ®CZ®...®CZ j ®I (24) 



JV/2-1 



All the unitary (|24|) implements are controlled phase gates CZ between qubits 2-3 and 
4-5 and so on. It connects the states p 2 together in a chain as shown in Fig. (J2j) . It is 
important to notice that although control-phase gates are entangling operations, they 
leave the state p* N in a separable form. 

This can easily be seen if we consider two separable 2-qubit states (j2J) joined by a 
control-phase operation. The total 4-qubit state can be written as pi = CZ 2 ^p\ <g> 
p 2 CZl 3 . This state is locally equivalent to p\ = |[|0000) (0000| + |0111)(0111| + 
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Figure 2. Closest separable state of N qubits. Each pair of qubits represents the 
closest separable state p^- 



|1011)(1011| + |1100)(1100|] where we have applied Hadamard operators on every second 
qubit / eg) H <g> / ® H . The transformed state pi is clearly separable. Since identical 
analysis can be applied to p* N we conclude that control-phase operations in the unitary 
(El]) do not introduce any entanglement and p* N remains separable. 

It is important to note that the form of the closest separable state is the same 
for cluster states in higher dimensions as well. The majority of the proof for closest 
separable states is independent of the dimension of our cluster state. The dimensionality 
becomes important only in Eq. (|16|) where we assume one-dimensional cluster state by 
using = CZ\ipk) 1^2) where the control-phase gate is applied between qubits k 

and k + 1 . For higher dimensions the number of control-phase gates in Eq. (fl6l) increases 
due to larger number of neighbouring qubits. However the logic of the proof remains 
unchanged. Therefore it is possible to find the closest separable state to a pure cluster 
state of any dimensionality. This is crucial because unlike linear cluster chains, two and 
higher dimensional clusters are universal resources for quantum computation [6]. 

3. Thermal entanglement 

Study of thermal entanglement is of great significance since any realistic scheme of 
implementing a quantum computer will operate at finite temperatures. In this section 
we look at entanglement in thermal cluster states. Specifically we are interested to see 
at what temperature the cluster state becomes separable. These considerations will lead 
us to a more general form for the closest separable state. 
Consider the cluster Hamiltonian 



where J is the coupling constant and M is the neighbourhood of site j. The terms in the 
sum above are just a particular case of stabilizer operators [6]. Hence Hamiltonian ( |25l ) 
is also referred to as the stabilizer Hamiltonian. The ground state of this Hamiltonian 
is the cluster state. Excited states are achieved by local cr| flips. The ground state has 
energy — NJ and is non-degenerate. The k th excited state has energy J(—N + 2k) with 
degeneracy k ^i k y [20] • All the states in the spectrum of the Hamiltonian are equally 
entangled since they are all some a z away from the ground state. Therefore our method 
of finding the closest separable state can also be used on these excited states. 



N 




(25) 
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The ground state of this Hamiltonian is highly entangled and the entanglement 
scales linearly [10] with the size of the cluster state as y. An interesting question to ask 
is whether this entanglement persists at finite temperatures and where is the critical 
point beyond which the thermal entangled state becomes separable. Thermal cluster 
state has the following form [12] 

1 N 

^)=^n[ /+tanh (w ( 26 ) 

3 

where Kj is the stabilizer operator at site j. For simplicity we will now consider a 
two-qubit thermal cluster state of the form 

a 2 (uj) = ^(I + uja x ® a z )(I + ujo z ® a x ) (27) 

where u = tanh(/3J). So in the low temperature limit as T — > we recover the pure 
state because to — > 1. As the temperature increases u — > and for T — > oo we obtain a 
maximally mixed state. To calculate the critical temperature at which the state becomes 
separable we use the Peres-Horodecki criterion [21] and [22]. To see when the state fails 
to be a positive definite we have to solve the following equation w\ + 2u c — 1 = 0. This 
equation has two solutions uj c = — 1 ± v2 but we will disregard the negative solution 
because it is not physical. Substituting this back into uj = tanh(/3J) finally gives us the 
critical temperature 

2 7 

T c = ^7= (28) 

fc B ln(V2-l) 

The entanglement for T < T c happens to be of useful distillable nature [23J. Therefore 
by using local operations and classical communication on multiple copies of the thermal 
cluster pure entanglement can be distilled. This critical temperature for distillable 
entanglement remains unchanged for higher dimensions and any system size iV as proved 
in [20]. For temperatures above the critical temperature, T > T c , the iV-qubit cluster 
state, where N > 2, is not fully separable but is bound entangled [21]. The critical 
temperature at which all entanglement vanishes is non-trivial to find due to the fact 
that properties of bound entanglement are still not fully understood. 

We have already calculated the closest separable state for a 2-qubit cluster. Now 
we will show that this state is not unique. In fact there are infinitely many states that 
satisfy condition (0). Since there is no bound entanglement for the case of two qubits 
we know that the thermal cluster state at the critical temperature T c is fully separable. 
It is natural to ask the question what the distance is between this state and a pure 
cluster state in terms of the relative entropy. So we want to calculate S^o^l |o"2(^c)) 
where a 2 (uj c ) = \{I + co c a x (g) cr z )(I + uj c & z ® cr x ) and tu c = v2 — 1. It can be quickly 
checked that S(a\ |cr 2 (ci; c )) = 1 and the gradient of relative entropy in a small region 
around cr 2 (k; c ) is |^(0, r) = 1 — 2Tr[ar] > as can be seen from Eq. fll4l) . Therefore the 
2-qubit thermal cluster at critical temperature is another good closest separable state 
to the pure cluster state. 



Entanglement in pure and thermal cluster states 



10 



(I<g,H)c 2 (u)(I<g>H^ = l - 



It is instructive to look at the thermal cluster state o 2 [uS) in the "mixed" basis and 
compare it to the pure cluster closest separable state p* 2 . Transforming 02(0;) we have 

f 1+u u(l + w)\ 

1-cj u{oj-1) 

u(u-l) 1-LU 

\u(l + u) 1 + u J 

From the form of the matrix it can be seen that the only way that <j 2 (uj) — > p 2 is when 
all the terms apart from 1+u vanish. So we arrive at the following four conditions that 
have to be satisfied simultaneously; lu(1 + uj) = 1 — u = u(u — 1) = and 1 + u = 2. 
It is straightforward to see that the only solution is u = 1 which corresponds to the 
pure cluster state case. Therefore as the temperature increases the thermal cluster state 
does not approach the pure closest separable state p 2 . This last result can be seen 
also from the fact that the relative entropy between these two states is non-zero, ie. 
S(p%\\(T2{wc)) 0- The relative entropy between two states vanishes if and 

only if A = B [26]. 

Taking a convex mixture of p 2 and o^^c) allows us to find a general form of the 
closest 2-qubit separable state 

p 2 = (l-\)p* 2 + \a 2 (u c ) (30) 

where A G [0, 1]. This can be verified by computing the gradient of the relative entropy 
HI 

Ox 



(29) 



in the usual way, §£(0, r) = 1 - 2Tr[ar] > for all A G [0, 1]. 



4. Thermal cluster states 



Now we are in a position to calculate the entanglement of thermal cluster states. We will 
use a similar approach as above to show that the form of the closest separable state is 
the same as in Eq. (l30|) . However this time the parameter A that determines the mixture 
of our two closest separable states from Eq. (1301) will not be completely free and its lower 
bound will depend on the temperature. We will show that as the temperature increases 
towards the critical temperature T — > T c the parameter A increases towards unity. 

We already have everything we need to compute the scaling of entanglement with 
temperature of a thermal 2-qubit cluster state given by Eq.( j27l) . The only thing that 
remains to determine is a suitable candidate for the closest separable state. The 
most obvious choice would be a 2 (u; c ) = Ml + uj c a x cg> cr z )(I + u c cr z ® a x ) where 
<jj c = \pl — 1. In fact this gives the correct behavior of entanglement E(a 2 (uj)) = 
Tr[a 2 (uj) log 02(0;) — cr 2 (u)) logcr 2 (a; c )] as can be seen in Fig. (j3J). The gradient of the 
relative entropy in the neighbourhood of the state ct 2 (uj c ) is non-neg ative ff(0,r) > 
as expected. 

Now that we know that for 2 qubits the cluster state at critical temperature is a 
closest separable state we can ask whether there exists a more general form of the state 
just like in Sec. (j3J). It turns out that this is in fact possible, albeit with some needed 
modifications. We have seen that as the temperature increases the thermal cluster state 
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Figure 3. Scaling of entanglement of 2 and 4 qubits with increasing temperature. 
When the thermal coefficient w = lwe recover the result for pure cluster state. At the 
critical point uj c = sfl — 1 the state is separable. 

does not approach the state p* N from Sec. ([2]). This can be also seen from the non- 
vanishing distance between the separable states <r N (u c ) and p* N , ie S(<jn(uj c )\\p* n ) ^ 0. 
As the temperature increases the distance between the cluster state (Tn(u)) and <tn(lj c ) 
approaches zero as expected. However this is not true for the distance between <Jn(uj) 
and p* N due to the fact that S(aiy(uj c )\\p* N ) ^ 0. 

This problem can be overcome by not allowing the parameter A from Eq. (l30l) to 
take any value from the interval [0,1]. Rather we require that the lower bound of the 
possible values of A increases with temperature. We call this bound A* and require 
A* = X*(u). 

As before we will first demonstrate the general principle at work on a two-qubit 
cluster state 02(0;). First we need to find the new parameter A* by minimising the 
distance between the thermal cluster state a2{uj) and the state (150]) P2 for a constant 
temperature. Solving J^S(a2(uj)\\p2) = for A gives the new A* 



(v / 2-2)(3 + cu) 

We can straight away confirm that the separable state given by this new parameter 
P2 = (1 — Ag)p2 + ^.oziyJc) is also a closest separable state by calculating the gradient 
of relative entropy at this state |£ (0, r) > 0. Therefore we can construct a new general 
closest separable state 

p 2 (u) = (l-\)p* + \a 2 (u c ) (32) 

where this time A G [AJj, 1]. Relative positions of all the above states in Hilbert space 
are illustrated in Fig.(TJ]). One thing that the figure does not capture is the fact that the 
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0"2 



E 



P2 = (1 - Aj)P2 + X^i^c) 



Figure 4. Section of Hilbert space of two qubits. The shaded area D represents 
the subspace of all separable states. E designates the set of all entangled states. As 
the temperature increases the set of all separable states becomes smaller. This set is 
represented by all the points found on the E-D boundary between the states <72 (w c ) 
and p2- 



distances S^Wp^) and S^c^l |cr 2 (c<; c )) are the same. This is because not all properties 
of Hilbert space can be pictured on a two-dimensional drawing. 

This method of finding closest separable states can be extended to general iV-qubit 
cluster states. The crucial part is to find the critical temperature at which the cluster 
becomes fully separable. Due to presence of bound entanglement for N > 3 this is a 
non-trivial task and at present it is unclear how it can be achieved. [23] calculates upper 
and lower bounds on the value of the critical temperature for cluster states of different 
dimensions. However all that this changes is that the separable state a(u c ) would be 
shifted to the left of the state a2(w c ) on the E-D boundary in Fig.(j3J) since the critical 
temperature for N > 3 is higher for than the one calculated here. 

5. Conclusion 

We have demonstrated a systematic way of constructing closest separable state to N- 
qubit cluster states. Our method reproduces known results for pure states and also 
allows us to quantify entanglement for mixed cluster states. 

For pure states the method relies on writing the density matrix of the state in 
the "mixed" basis and setting any off-diagonal elements to zero. Operationally we 
can construct this state from N/2 copies of 2-qubit closest separable state and joining 
them together by applying control-phase gates as demonstrated in Sec.(j2D- Our method 
applies to cluster states of all dimensions. This is particularly useful since cluster states 
of dimension d > 2 are universal resource state for quantum computation. 

It also turns out that state p* N is not the only closest separable state. Thermal 
cluster state at critical temperature (Jn{oj c ) is also another closest separable state to 
pure cluster state. Therefore any convex mixture of these two states also minimizes 
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the relative entropy. However at the moment it is difficult to determine the precise 
temperature at which the cluster becomes fully separable for cases of more than 2 
qubits. This is due to the presence of bound entanglement. 

Mixed states require an even more careful approach. When quantifying 
entanglement in thermal cluster states a convex mixture of p* N and a^iojc) does not 
minimize the relative entropy anymore. Parameter A that determines the mixture of 
these two states does not take any value from interval [0, 1]. Instead as the temperature 
increases the lower bound of this interval increases as well. So the range of values the 
parameter A can take is constricted to A e [A*, 1]. The new lower bound A* is a function 
of temperature and as temperature approaches the critical value A* — > 1. In other words 
both the thermal cluster and the closest separable state p approach the same state. 
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